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Weak corrections and high Et jets at Tevatron 
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We calculate one-loop (purely) Weak (W) corrections of O{oi^oofi) to the partonic cross section of 
two jets at Tevatron and prove that they can be larger than the tree-level ©(qsOew) and C1 (q|;w) 
Electro-Weak (EW) ones. At high transverse energy of the jets, all such corrections may lead to 
detectable effects of, e.g., —10% or so, with respect to the leading-order (LO) QCD term of 0{ag), 
for the highest value so far probed by Run 2, depending on the factorisation/renormalisation scale. 
Besides, they increase significantly with jet transverse energy. Hence, our results show that EW 
corrections may be needed to fit the Standard Model (SM) to present and future Tevatron jet data. 


As the overall energy of hard scattering processes in¬ 
creases one should expect a relatively large impact of 
perturbative EW corrections, as compared to the QCD 
ones. This can easily be understood (see [1] - [2] and ref¬ 
erences therein for reviews) in terms of the so-called Su- 
dakov (leading) logarithms of the form a^w log^{s/M^) 
(hereafter, a-w = Q;em/ sin^ 0w, with oem the Electro- 
Magnetic (EM) coupling constant and the weak mix¬ 
ing angle, whereas >/! is the parton-level centre-of-mass 
energy), which appear in the presence of higher order 
weak corrections when the initial state carries a definite 
non-Abelian flavour and which, unlike QCD, do not can¬ 
cel between virtual and real emission of W bosons [3]. 

Furthermore, one should recall that real weak bosons 
are unstable and decay into high transverse momentum 
leptons and/or jets, which are normally captured by the 
detectors. In the definition of a hadronic cross section, 
one may then remove events with such additional parti¬ 
cles. Hence, for typical experimental resolutions, softly 
and collinearly emitted weak bosons need not be included 
in the definition of the production cross section and one 
can restrict oneself to the calculation of weak effects orig¬ 
inating from virtual corrections only. In fact, leading 
(and all subleading) virtual weak corrections are finite 
(unlike QCD, where infrared divergences mean that vir¬ 
tual corrections must be considered in conjunction with 
gluon bremsstrahlung), as the mass of the weak gauge 
boson provides a physical cut-off for the otherwise diver¬ 
gent infrared behaviour. Under these circumstances, the 
(virtual) exchange of Z bosons also generates double log¬ 
arithmic corrections, aw {s/M"^). Moreover, in some 
simpler cases, the genuinely weak contributions can be 
isolated in a gauge-invariant manner from purely EM ef¬ 
fects and the latter may or may not be included in the 
calculation, depending on the observable being studied. 

The leading, double-logarithmic, angular-independent 
weak logarithmic corrections are universal, i.e., they de¬ 
pend only on the identities of the external particles. 
In some instances, however, large cancellations between 
angular-independent and angular-dependent corrections 
[4] (see also [5] for two-loop results) and between leading 
and subleading terms [6] have been found at TeV ener¬ 
gies. Moreover, some other considerations are in order 
in the specihc hadronic context. Firstly, one should re¬ 


call that hadron-hadron scattering events involve valence 
(or sea) partons of opposite isospin in the same process, 
but since the PDFs are not singlets of flavour only par¬ 
tial cancellations among initial state large logarithms will 
occur [3]. Secondly, several crossing symmetries among 
the involved partonic subprocesses can also easily lead to 
more cancellations. 

Because all this, it becomes of crucial importance to 
study the full set of fixed order weak corrections, in view 
of establishing the relative size of the different contribu¬ 
tions at the energies which can be probed at TeV scale 
hadronic machines. Several results already exist, e.g., in 
the SM, for: EW gauge boson production in single mode 
[4], [7] as well as in pairs [8]; bb [9] and tt [10] - [15] 
production; Higgs processes [16]. (See [17] for a review.) 

It is the aim of our paper to report on the computation 
of the full one-loop weak effects* entering all possible ‘2 
parton ^ 2 parton’ scatterings, through the perturba¬ 
tive order agaw- (See Ref. [18] for tree-level asaEw 
interference effects - hereafter, oew exemplifies the fact 
that both EM and W effects are included at the given 
order). We will ignore altogether the contributions of 
tree-level agaw terms involving the radiation of W and 
Z bosons. Therefore, apart from gg —*■ gg, qq' QQ', 
qq' — *■ QQ’ and qq' — *■ QQ' (which are not subject to 
order agaw corrections), there are in total fifteen sub¬ 
processes to consider, 

99 qq^99, (i) 

q9 qq, qq qg, (2) 

qq qq, qq qq, (3) 

qQ —!■ qQ (same or different generation), (4) 

qQ — !■ qQ (same or different generation), (5) 

qq qq, (6) 

qq QQ (same or different generation), (7) 

qQ qQ (same or different generation), (8) 


*We neglect considering here purely EM effects (as well as 
interferences between these and the weak ones), as they can 
be isolated in a gauge invariant fashion and since they are not 
associated with logarithmic enhancements either (like QCD). 
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with ) and ^ referring to quarks of different flavours, 
limited to u-, d-, s-, c- and 6-type (all massless). While 
the first four processes (with external gluons) were al¬ 
ready computed in Ref. [19], the eleven four-quark pro¬ 
cesses are new to this study (see Ref. [20] for RHIC and 
LHC results). Besides, unlike the channels with exter¬ 
nal gluons, those with four-quarks must include virtual 
gluon corrections to tree-level interferences between weak 
and strong interactions and therefore can be infrared di¬ 
vergent, which means that gluon bremsstrahlung effects 
must be evaluated to obtain a finite cross section at the 
given order. In addition, for completeness, we have in¬ 
cluded the non-divergent subprocesses of (anti-)quark- 
gluon scattering into three coloured fermions. 

Our studies are of particular relevance in the context of 
the Tevatron collider, where an excess was initially found 
by CDF (but not DO) at high transverse energy in inclu¬ 
sive jet data from Run 1 [21], with respect to the next- 
to-LO (NLO) QCD predictions [22] - [24]. While sev¬ 
eral speculations were made about the possible sources 
of such excess from physics beyond the SM, it was even¬ 
tually pointed out that a modification of the gluon PDFs 
at medium/large Bjorken x can apparently reconcile the¬ 
ory and data within current systematics: see, e.g., [25]. 
(For a different explanation, see [26].) In fact, notice 
that with the most recent PDFs (e.g., CTEQ6.1M [27]), 
also preliminary Run 2 data seem to be (barely) consis¬ 
tent with NLO QCD, see [28] for CDF. (Results from 
DO have a larger systematic uncertainty, which tends to 
encompass the theory predictions [28].) 

Over a hundred one-loop and tree-level diagrams are 
involved in the computation of processes (l)-(8) and is 
thus of paramount importance to perform careful checks. 
In this respect, we should mention that our expressions 
have been calculated independently by at least two of us 
using FORM [29] and that some results have also been 
reproduced by another program based on FeynCalc [30]. 

As already mentioned, infrared divergences occur when 
the virtual or real (bremsstrahlung) gluon is either soft 
or collinear with the emitting parton and these have 
been dealt with by using the formalism of Ref. [31], 
whereby corresponding dipole terms are subtracted from 
the bremsstrahlung contributions in order to render the 
phase space integral free of infrared divergences. The 
integration over the gluon phase space of these dipole 
terms was performed analytically in d-dimensions, yield¬ 
ing pole terms which cancelled explicitly against the pole 
terms of the virtual graphs. There remains a divergence 
from the initial state collinear configuration, which is ab¬ 
sorbed into the scale dependence of the PDFs and must 
be matched to the scale at which these PDFs are ex¬ 
tracted. Through the order at which we are working, it 
is sufficient to take the LO evolution of the PDFs (and 
thus the one-loop running of as)- 

Some of the diagrams also contain ultraviolet diver¬ 
gences. These have been subtracted using the ‘modi¬ 


fied’ Dimensional Reduction (DR) scheme at the scale 
/i = Mz- The use of DR, as opposed to the more usual 
‘modified’ Minimal Subtraction (MS) scheme, is forced 
upon us by the fact that the W- and Z-bosons contain 
axial couplings which cannot be consistently treated in 
ordinary dimensional regularisation. Although not essen¬ 
tial, we find it convenient to work with helicity matrix el¬ 
ements extracted using properties of Dirac matrices valid 
in four dimensions. Thus the values taken for as refer 
to the DR scheme whereas the EM coupling, aEM, has 
been taken to be 1/128 at the above subtraction point. 
(The numerical difference between these two schemes is 
negligible for as though.) 

For the top mass and width, entering some of the loop 
diagrams with external 6-quarks, we have taken rrit = 
175 GeV and Ft = 1.55 GeV, respectively. The Z mass 
used was Mz = 91.19 GeV and was related to the W 
mass, Mw, via the SM formula Mw = Mz cosOw, where 
siv? 6w = 0.232. (Gorresponding widths were Tz = 2.5 
GeV and Fw = 2.08 GeV.) 


jet-production (0.1 < [rj] < 0.7) 
6 = (NL0-L0)/L0 
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FIG. 1. The effects of the G(Qgaw) corrections rel¬ 
ative to the LO SM results for Run l(Run 2) using 
CTEQ3L(CTEQ6L1) as PDFs. 

Fig. 1 shows the effects of our one-loop corrections 
to the LO results for jet production, the latter be¬ 
ing defined as including all possible terms of order a|, 
Q!sc«EW and (hereafter LO SM). (The spike at 

Et ~ Mwl‘^,Mzl‘2. is a threshold effect in the loop di¬ 
agrams.) Notice that in our treatment we identify the 
jets with the partons from which they originate and we 
adopt here the cut 0.1 < \ri\ < 0.7 in pseudorapidity to 
mimic the GDF detector coverage and the standard jet 
cone requirement Ai? > 0.7 to emulate the jet data selec¬ 
tion (although we eventually sum the two- and three-jet 
contributions). Furthermore, as factorisation and renor¬ 
malisation scale we use ^ = fip = Et/2 - a choice 

leading to the best convergence of both NLO [24] and 
resummed [32] QGD predictions - (where Ep is the jet 
transverse energy) while we adopt GTEQ3L as PDFs [27] 


2 




for Run 1, a set defined prior to the re-arrangement of the 
gluon. With respect to the LO SM rates, the 0(Q;gaw) 
corrections are not large despite growing steadily with 
Et ■ For Et values in the vicinity of 420 GeV, the highest 
point of Run 1 and also the location of the apparent CDF 
excess, they amount to —1.5%. This effect is not com¬ 
petitive with the positive NLO QCD corrections through 
0{ag): see, e.g.. Fig. 1 of [24]. In the same figure, we 
have also shown the ©(agcrw) corrections at Run 2 for 
the same fjL and the choice CTEQ6L1 of PDFs (one of the 
newest sets incorporating the above mentioned gluon re- 
parameterisation). Here, we have also increased the Et 
values probed, as the larger collider energy has already 
allowed to collect data some 150 GeV beyond the Run 1 
reach. We see that at the higher energy the 0(Q;gaw) 
corrections are substantially similar in size and shape to 
the lower energy case, so that they stretch to —2% or so 
near the current kinematic limit (550 GeV or so). (Cross¬ 
ing points between the two curves are induced by the dif¬ 
ferent PDF choice as well as the different numerical value 
of ^ at the two energies.) 


jct-production ([tj] < 2.5) 
6 = (NLO-LO)/LO 



0 200 400 600 800 

Et (GeV) 

FIG. 2. The effects of the (!l(a|aw) corrections relative to 
the LO SM results for Run 2 in the presence of two sets of 
up-to-date PDFs (CTEQ6L and CTEQ6L1). 

Fig. 2 extends the Et interval to 850 GeV and the 
pseudorapidity cover to jryj < 2.5 (our new default from 
now on, for same Ai?), while still adopting /i = Et/‘2. 
as factorisation/renormalisation scale. Including the for¬ 
ward/backward detector region reduces minimally the ef¬ 
fects of the ©(agaw) corrections while their shape re¬ 
mains unchanged. Their maximum is about —5% at the 
upper end of the interval considered. Furthermore, their 
dependence on the choice of PDFs is also very small, as 
we have verified by running CTEQ6L1 vs. CTEQ6L [27]. 

Notice however that, if one defines the corrections with 
respect to only the 0{a^) contribution (hereafter, LO 
QGD), the effects of the sum of all non-QCD terms, i.e., 
those of order ®ew Q^s'^w (hereafter LO SM 

-I- NLO W), become significantly larger. Fig. 3 makes 


this point clear. At Et = 850 GeV or so, the upper 
kinematic limit of the collider, one would see a combined 
effect of about —14%, most of which are indeed due to the 
O(agaw) terms new to this study (NLO W). In practice 
though, such jet transverse energies are unreachable even 
for optimistic luminosity. For the current Run 2 highest 
Et point, 550 GeV, the effects of the LO SM -I- NLO 
W corrections amount to —8% of the LO QGD term. 
Clearly, it is of paramount importance to establish which 
terms are included in Monte Carlo (MC) programs used 
to interpolate the data. In general, it is clear from Fig. 3 
that the corrections due to the one-loop graphs play a role 
at least as relevant as those due to tree-level effects and, 
importantly, at Tevatron, they act in the same direction, 
namely, a reduction of the differential QCD rates. 


jet-production {[tj] < 2.5) 
5 = (NL0-L0)/L0 
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FIG. 3. The effects of the ©(osoew -I-Obw) (LO SM) and 
the latter plus 0{a%aw) (LO SM -|- NLO W) corrections 
relative to the LO QCD of 0{ag) results for Run 2 in the 
presence of up-to-date PDFs (CTEQ6L1). 

In fact, another subtlety should be borne in mind as 
far as EW corrections are concerned. We have so far 
adopted /i = Et/‘2. for the factorisation/renormalisation 
scale. This seems in fact to be the preferred choice while 
comparing Tevatron data against NLO QCD predictions 
through (!I(ag). A discussion of the dependence of the 
QCD corrections on /r is found in Refs. [23] - [24] and 
the above mentioned choice is motivated by the stability 
of the higher order QCD results in the region /r Et/‘2- 
In fact, recall that any dependence on ^ arises because 
of the truncation of the perturbative expansion at some 
fixed order and it is therefore a measure of the missing 
higher order terms. As /r would not appear if these were 
known through all orders, it is customary to vary the 
factorisation/renormalisation scale in order to estimate 
the residual theoretical error. We have done so in Eig. 4 
for, e.g., Et = 100 and 550 GeV, at Run 2 energy with 
CTEQ6L1 as PDEs. The fact that the C>(agaw) curves 
do not display local maxima, unlike the 0(ag) results 
(Fig. 2 of [23]), does intimate that one scale choice is 


3 










not more appropriate than another (irrespective of the 
jet transverse energy probed and the size of the EW cor¬ 
rections). Thus, there is no firm reason to adopt Et!^ as 
factorisation/renormalisation scale here. If a higher value 
is chosen at 550 GeV, e.g., ^ = Et, the LO SM -I- NLO 
W corrections grow of a further percent, to —9%, while 
for /i = 2Et they become —10%. This trend is manifest 
over the entire Et range of relevance at Tevatron. 


jet-production (|7j| < 2.5) 

V s = 2000 GeV 6 = (NL0-L0)/L0 /x = ^ (GeV) 



FIG. 4. The effects of the ©(osobw + olw + oiow) (LO 
SM -I- NLO W) corrections relative to the LO QCD results as 
function of p for Run 2 in the presence of up-to-date PDFs 
(CTEQ6L1) for two choices of jet transverse energy. 

In summary, at Tevatron, EW effects in general and 
0(a|aw) one-loop terms in particular are important 
contributions to the inclusive jet cross section at large 
transverse energy. A careful re-analysis of actual jet data, 
which was beyond the intention of this paper, may be 
needed in view of the increasing luminosity of the Fer- 
milab collider. Particular care should be devoted to the 
treatment of real W and Z production and decay in the 
definition of the inclusive jet data sample, as this will 
determine whether tree-level W and Z bremsstrahlung 
effects have to be included in the theoretical predictions 
through (P(a|aw), which might counterbalance the neg¬ 
ative effects due to one-loop W and Z virtual exchange. 
In closing, we should mention that the calculation of the 
aforementioned EM effects is in progress. 
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